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ABSTRACT 

The  purpose  of  this  thesis  is  to  present  solution  tech¬ 
niques  for  provisioning  problems  arising  in  the  Navy’s 
wholesale  purchase  of  replacement  parts.  The  objective  is 
to  minimize  the  Mean  Supply  Response  Time  (MSRT)  subject  to 
a  budget  constraint.  The  problem  can  be  formulated  as  a 
Dynamic  Program  (DP) ,  however,  it  is  too  large  and  complex 
for  a  standard  recursive  dynamic  approach.  Consequently,  a 
variation  of  the  normal  DP  approach  was  developed  that 
significantly  reduces  the  required  computations.  An 
existing  DP  computer  program  was  modified  to  implement  this 
DP  variation.  The  result  is  a  usable  approach  considering 
speed  and  ease  of  manipulation. 
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I.  HISTORY  AND  IN TROD OCT ION 


The  U.  S.  Navy  is  a  hardware  oriented  organization 
requiring  a  wide  variety  of  spare  and  replacement  parts. 
The  Naval  Supply  Systems  Command  (NA  VS  UP)  is  the  organiza¬ 
tion  within  the  Navy  responsible  for  maintaining  a  suffi¬ 
cient  stock  of  parts  and  replenishing  this  stock  when 
quantities  run  low.  NAVSOP  is  keenly  interested  in  the 
guality  of  the  models  it  uses  to  provide  support  for  the 
Navy's  weapons  systems.  Faculty  of  the  Naval  Postgraduate 
School  were  asked  to  examine  the  existing  wholesale  and 
retail  provisioning  and  replenishment  models  and  to  provide 
comments  on  their  correctness  and  suggestions  for  possible 
improvements.  In  respoa se#  a  report  by  Richards  and 
McMasters  [Ref.  1]  addresses  the  problems  and  various  solu¬ 
tion  techniques.  This  thesis  examines  a  specific  solution 
technique  that  may  be  applicable  to  one  of  the  problems. 
The  following  is  a  brief  synopsis  of  the  existing  models  and 
the  environment  in  which  they  operate. 

Within  NA  VS  UP  there  are  two  centers  that  use  provi¬ 
sioning  models,  1)  Ships  Parts  Control  Center  (SPCC)  and  2) 
Aviation  Supply  Office  (ASO) .  Presently  the  Inventory 
Control  Points  (ICPs)  within  the  centers  have  two  retail 
provisioning  models  which  generate,  1)  COSALS  (Coordinated 
Shipboard  Allowance  Lists)  and  2)  AVCALS  (Aviation 
Coordinated  Allowance  Lists)  .  These  models  supply  suffi¬ 
cient  parts  support  for  a  short  term  ship  deployment  and  try 
to  provide  protection  from  replacement  part  depletion  until 
the  wholesale  reple rishmen t  models  can  buy  or  repair  the 
needed  spare  parts. 

Weapon  system  procurement  requires  that  provisioning 
support  be  available  by  a  set  preliminary  operational 


capability  (POC)  date.  The  initial  support  requirements  are 
determined  by  a  wholesale  provisioning  model.  The  support 
must  be  comprehensive  enough  to  supply  a  new  weapon  system 
until  the  replenishment  buy  can  be  made  and  delivered.  It 
should  be  noted  that  the  provisioning  lead  time  is  usually 
two  to  two  and  one  half  years.  Variability  in  production 
and  installation  causes  the  quality  of  the  forecasted  demand 
to  be  poor.  Hence,  the  provisioning  buy  may  result  in  large 
quantities  of  excess  parts  in  the  supply  system,  or  there 
may  be  many  stockouts  before  the  replenishment  buy  becomes 
available.  The  Assistant  Secretary  of  Defense 
(Installations  and  Logistics)  has  provided  guidance  [Ref.  2] 
to  help  prevent  large  excesses.  From  this  instruction  a 
third  model  (COSDIF)  was  developed  to  conservatively  esti¬ 
mate  the  quantities  of  new  parts  required  for  wholesale 
provisioning.  For  existing  equipment  no  additional  stock  is 
procured  for  the  new  weapons  system.  If  the  increase  in 
demand  from  the  new  system  is  large  enough,  the  inventory 
manager  can  anticipate  the  increase  and  start  a  replenish¬ 
ment  buy  in  sufficient  time  to  meet  the  need.  Otherwise, 
inventory  optimization  models  may  be  used  to  control  the 
stock.  Wholesale  replenishment  for  a  new  system  is  initi¬ 
ated  at  different  times  by  the  two  centers  within  NAVSUP. 
For  either  center  the  reorder  point  for  replacement  parts  is 
subject  to  change  as  a  function  of  the  actual  demand  the 
weapon  system  places  on  the  supply  system.  The  reorder 
point  for  an  existing  system  occurs  at  part  depletion  levels 
set  by  NAVSOP. 

Budgeting  constraints  for  support  of  existing  systems 
are  partially  driven  by  the  peculiarities  of  the  Department 
of  Defense  Planning,  Programming  and  Budgeting  System 
(PPBS).  Those  purchases  for  new  weapon  systems  are  funded 
as  part  of  the  procurement  process.  For  existing  systems, 
each  Service  requests  from  Congress  sufficient  funding  to 
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of  costs.  Consequently,  a  variation  of  the  normal  D? 
approach  was  developed  that  significantly  reduces  the 
required  computations.  This  variation  shall  be  referred  to 
as  'the  funnel'  because  of  its  funneling  characteristics  for 
computation  of  the  resource  variable  bounds.  Further  expla¬ 
nation  of  its  construction  and  procedure  may  be  found  in 
Chap  2S3.  Modifications  of  an  existing  DP  computer  program 
described  in  [  Eef.  3]  to  operate  as  the  DP  variation  makes 
the  'funnel'  method  a  usable  approach  considering  speed  and 
ease  of  manipulation.  Additionally,  tne  DP  variation  has 
the  advantage  of  providing  an  optimal  integer  solution.  In 
this  report  'the  funnel'  will  be  compared  to  a  full  DP 
procedure  and  a  computeri red  marginal  analysis  approach 
which  operates  very  efficiently  but  which  may  not  generate 


II.  SOLOTION  TECHNIQOES 


This  chapter  will  discuss  three  methods  foe  solving  the 
provisioning  problem  described  in  Chapter  1,  1)  marginal 

analysis,  2)  pure  dynamic  programming,  and  3)  a  modified 
dynamic  programming  procedure.  The  first  technigue  provides 
a  guick  solution  that  is  a  good  estimate  but  may  not  be 
optimal.  The  second  technique  provides  a  global  optimal 
solution  but  reguires  a  large  amount  of  computer  time  for 
any  moderate  sized  problem.  The  third  technigue  provides  a 
local,  and  possibly  a  global,  optimal  solution  and  reguires 
significantly  less  time  than  option  two. 

The  following  definitions  will  be  used  in  the  formula¬ 
tion  and  discussions  in  this  chapter.  Capital  letters 
represent  the  variable  while  lower  case  letters  indicate  a 
subscript.  In  the  expression  Xi-1  the  quantity  i- 1  is  the 
subscript  although  it  appears  on  the  same  line  as  the  symbol 
X. 

N  =  the  total  number  of  items  considered  for 

provisioning 

3  =  pre-assigned  maximum  budget  amount 

Ci  =  the  unit  cost  of  item  i 

Ei  =  the  essentiality  coda  for  item  type  i 

Hi  =  the  fixed  time  required  to  satisfy  a 
requisition  for  item  i  if  stock  is 
available 

Li  =  the  demand  rate  for  item  i  (lambda) 

Si  =  the  decision  variable  for  the  number 

of  items  of  type  i 

Ti  =  the  procurement  lead  time  for  item  i 
P'i(Si)=  the  probability  of  Si  demands  for  item  i 
during  the  provisioning  interval 
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Pi  (Si)  =  the  protability  of  Si  or  fewer  demands 
for  item  i  daring  the  provisioning 
interval 

Zi  (Si)  =  the  performance  measure  for  item  i  when 
Si  units  are  stocked. 

The  objective  of  this  problem  is  to  minimize  the  Mean 
Supply  Eesponse  Time  (MSRT)  as  a  function  of  Si.  M SRT  is  a 
widely  used  measure  cf  supply  performance  because  of  its 
role  in  determining  availability  and  because  it  is  an  indi¬ 
cator  of  the  success  of  a  supply  system  in  meeting  response 
time  goals.  It  considers  two  major  factors  in  computation, 

1)  the  likelihood  of  satisfying  demands  from  stock  on 
hand  and 

2)  the  length  of  the  delay  in  satisfying  demands 
when  the  system  runs  out  of  stock. 

It  may  be  represented  by  the  following  expression, 

MSRT=min  £Ei*Li*Ti*Zi  (Si)  /£Ei*Li*Ti  (2.1) 

1  =  1  ,  •  .  .  ,  N  , 

where  the  likelihood  adjustments  for  conditions  one  and  two 
are  factors  in  the  Zi(Si)  function.  This  function  repre¬ 
sents  the  performance  measure  of  item  i  when  Si  items  are 
considered  and  is  defined  by  the  following  equation  for  the 
case  where  demands  during  the  provisioning  leadtime  are 
assumed  to  be  Poisson  distributed: 

Zi  (Si)  =Ki+  ( (1-Pi  (Si))*(Li*Ti-2Si*Si*(Si*1)  (2.2) 

/Li*Ti)  *P'  i  (Si)  (Li*Ti-Si)  )  /2Li)  . 

The  constraint  for  this  problem  is  the  budget  or  the 
resource  in  dollars.  It  is  linear  and  may  be  represented  by 
the  following  inequality. 
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£  Ci*si<=E 


i=  1 


t  •  •  •  / 


(2.3) 


N. 


The  following  subsections  present  solution  techniques  which 
are  viable  procedures  for  solving  the  above  formulated 
provisioning  problem. 

A.  MARGINAL  ANALYSIS 

The  marginal  analysis  method  is  a  quick,  simple  and 
effective  approach  for  generating  a  solution  to  the  previ¬ 
ously  defined  provisioning  problem.  For  the  problem  consid¬ 
ered  here,  the  marginal  analysis  approach  generates 
solutions  which  may  not  be  optimal.  However,  the  solutions 
generated  by  the  procedure  do  possess  some  very  desirable 
properties  which  are  stated  below. 

Let  Zi(Si)  be  the  essentiality  weighted  MSRT  for  item  i 
when  Si  units  are  stocked.  Let  &Zi  (Si)  =Zi  (Si- 1 )  - Zi  (Si)  be 
the  reduction  in  essentiality  weighted  MSRT  achieved  by 
increasing  the  stockage  level  for  item  i  from  Si-1  units  to 
Si  units.  Let  S' =  (S1,S2,.  . . , SN)  be  the  vector  of  alloca¬ 
tions  for  tin  N  items:  i.e.  Si  is  the  number  of  units  allo¬ 
cated  to  item  i.  The  marginal  approach  begins  with  a  budget 
of  B  dollars  and  an  allocation  vector  S' *  (0,0, . .  . ,  0)  .  It 
computes  £Zi(1) /Cl  for  i=1,2,...,N  and  selects  that  item  for 
which  this  marginal  benefit  to  cost  ratio  is  greatest. 
Without  loss  of  generality  ,  assume  that  the  item  selected 
is  item  1.  This  requires  an  expenditure  of  Cl  dollars  and 
results  in  the  allocation  vector  S  ’  =  ( 1 , 0, 0  ,. .  .  ,  0)  .  The 
procedure  next  computes  £Z1(2)/C1,  the  marginal  benefit  to 
cost  ratio  achieved  by  increasing  the  stockage  for  item  1  by 
one  unit  from  a  level  of  one  to  two  units.  The  maximum 
value  is  then  selected  from  among  the  ratios  in  2.4. 

An  additional  unit  of  the  item,  say  item  j,  for  which  the 
maximum  occurs  is  selected  and  the  remaining  budget  is 


(2.4) 


Az  1 ( 2 )  ,  Az 2 ( 1 )  ,  Az3(1) 

Cl  C  2  C  3 


,  AZN  (  1  ) 
CN 


decremented  by  C j  dollars.  The  process  continues  in  this 
manner  allocating  a  single  unit  iteratively  until  the  budget 
is  expended  or  no  additional  units  can  be  purchased.  At 
each  step,  the  ratios  compared  are 


Az  1  (S  1  +  1 )  ,  Az2  (S2  +  1 )  , 

C  1  C  2 


.  ,  AZN  (  SN+1) 
CM 


(2.5) 


where  the  current  allocation  vector  is  S'  =  (S  1 ,  S2 , . .  .SN)  . 
Since  only  one  of  the  ratios  changes  at  each  step  of  the 
process,  the  computations  required  are  very  few. 

It  can  be  shown  [Ref.  4],  that  for  the  problem  consid¬ 
ered  here  the  solution  is  undominated.  That  is,  if  B(S')  is 
the  budget  required  tc  fund  the  allocation  vector  S',  and  if 
S”  is  any  other  allocation  vector  such  that  3  (S'  • )  <=3  (S' ) 
then  the  overall  essentiality  weighted  HSRT  for  S'  is 
smaller  than  that  for  S**.  If  it  should  happen  that  B(S')  = 
B,  the  original  budget,  for  some  step  of  the  marginal  anal¬ 
ysis  procedure,  then  the  solution  S'  is  optimal. 

In  actual  practice,  what  usually  happens  in  the  marginal 
analysis  procedure  is  that  at  some  step  purchase  of  the  item 
with  the  largest  marginal  benefit  to  cost  ratio  requires  an 
expenditure  of  more  money  than  that  which  remains.  The 
allocation  vector  at  that  point  is  not  feasible  and  various 
heuristics  are  usually  applied  to  continue  the  process.  The 
most  popular  heuristic  is  to  backtrack  to  the  previous 
feasible  solution  and  select  that  item  which  has  the  maximum 
benefit  to  cost  ratio  among  those  items  with  unit  costs  no 
larger  than  the  budget  remaining.  On  applying  such  heuris¬ 
tics,  one  can  not  claim  that  the  final  solution  so  obtained 
is  undominated.  However,  one  can  make  use  of  the  undomi¬ 
nated  property  to  obtain  a  bound  on  the  optimal  value  of  the 
objective  function. 
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Let  Z(SS)  be  the  objective  function  value  at  the  occa¬ 
sion  at  which  the  marginal  analysis  procedure  first  gener¬ 
ates  a  solution  for  which  the  budget  required  is  equal  to  or 
greater  than  B.  If  equality  exists,  SS  is  optimal.  Suppose 
E(SS)=BS>B  and  let  S*  (X)  be  the  optimal  solution  for  a 
budget  of  X  dollars.  Then,  since  SS  is  undominated, 
Z(SS)  <=Z  (S*  (BS)  )  <=Z  (S*(B) )  ,  which  is  the  value  of  the  objec¬ 
tive  function  at  the  optimal  solution.  Thus  Z(SS)  is  a 
lower  bound  on  the  optimal  value  of  the  objective  function. 
Cne  can  thus  assess  the  potential  benefit  of  a  more  exten¬ 
sive  search  for  the  optimal  solution  by  comparing  the  actual 
performance  obtained  at  any  stage  with  the  lower  bound 
discussed  above. 

A  computer  program  was  written  by  Richards  to  implement 
this  method.  A  sample  of  this  Marginal  Analysis  program 
output  in  table  I  shows  the  iteration  number,  which  items 
are  purchased,  the  numbers  of  items  purchased,  the  benefit 
to  cost  ratios  and  the  budget  remaining.  The  table  shows 
the  sequence  in  which  the  items  are  purchased.  Note  the 
decreasing  benefit  to  cost  ratio.  The  program  selects  item 
number  one  five  times  in  a  row  before  the  first  unit  of  item 
two  provides  a  tetter  marginal  return.  The  last  attempt  to 
purchase  chooses  an  item  that  drives  the  remaining  budget 
negative.  The  program  automatically  re-evaluates  the 
marginal  returns  and  then  chooses  the  next  largest  benefit 
to  cost  ratio  which  in  this  case  drives  the  remaining  budget 
to  zero.  As  previously  discussed  this  solution  may  not  be 
optimal. 

The  Marginal  Analysis  program  has  several  other  capabil¬ 
ities  based  on  the  same  approach.  There  is  a  job  selection 
menu  that  provides  a  choice  of  1)  spares  determination,  2) 
budget  determination  or  3)  kit  evaluation.  The  first  option 
selects  the  quantities  to  be  purchased  within  a  given 
budget.  The  second  option  solves  a  related  problem  by 
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TABLE  I 

Sample  Output  from  the  Marginal  Analysis  Program 


*********************************************** 
*  METHOD:  FIXED  EUDGET  =  205-00  * 

*********************************************** 

*********************************************** 


STEP 

ITEM 

SI 

RATIO 

BR 

1 

1 

1 

0.480808496 

200.00 

2 

1 

2 

0.365659833 

195.00 

3 

1 

3 

0.260618091 

190.00 

4 

1 

4 

0. 172421157 

185.00 

5 

1 

5 

0.  105280340 

180.00 

6 

2 

1 

0.0735758543 

170.00 

7 

3 

1 

0.0600002967 

155.00 

8 

1 

6 

0.0591956675 

150.00 

9 

3 

2 

0.0533369593 

135.00 

10 

3 

3 

0.0466887541 

120. CO 

11 

3 

4 

0. 04C0909968 

105. 00 

12 

3 

5 

0.0336193480 

90.00 

13 

1 

7 

0.0306577347 

85.00 

14 

3 

6 

0.0273999237 

70.00 

15 

3 

7 

0.0216008648 

55.00 

16 

2 

2 

0-0207276568 

45.00 

17 

3 

8 

0.0164023377 

30-00 

18 

1 

8 

0.0146531016 

25.00 

19 

3 

9 

0.0119544677 

10.00 

20 

3 

10 

0.0083406679 

-5.00 

20 

1 

9 

0.0064818673 

5.00 

21 

1 

10 

0.0026624922 

0.0 

a***************************************** 

MINIMIZE 

MSRT 

EXCESS:  0.0 

I 

SI 

MSRT 

1 

10 

0.  190 

2 

2 

10.396 

3 

9 

12- 176 

OVEfi ALL  5. 385 

REALIZED  PERFORMANCE  =  5.3852 

BOUNDS  (  6.0033,  4.3120) 


computing  the  required  budget  to  meet  a  specified 
constraint,  in  this  case  MSRI.  The  third  option  evaluates 
the  available  performance  for  a  given  decision  set.  The 
program  will  conduct  the  optimization  problem  with  a 
laSrange  Multiplier  or  marginal  analysis.  In  either  case 
the  answer  is  the  same.  The  speed  and  ease  of  manipulation 
of  this  procedure  are  desirable;  however,  to  be  guaranteed 
an  optimal  solution  a  means  of  verification  must  be  found. 
That  is  the  goal  of  this  thesis. 

B.  STANDARD  DYNAMIC  PROGRAMMING 

The  standard  dynamic  programming  approach  can  be  applied 
to  this  type  of  provisioning  problem.  The  following  sample 
problem  illustrates  how  the  concept  may  be  applied.  A 
single  stage  problem  corresponding  to  a  single  item  may  be 
characterized  as  in  fig  2.1. 


Y  =  T ( X , S  ) 


Figure  2.1  Example  of  a  One-Stage  Decision  Problem. 

In  figure  2.1, 

X  is  the  input  state  variable  which  is  the  resource 
available  for  use  at  this  stage. 

S  is  the  decision  variable  that  represents  the  decision 
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within  the  stage. 


E  is  the  stage  return  function  that  produces  a  return 
based  on  the  input  X  and  the  decision  S. 

X  is  the  output  state  which  is  a  function  of  X  and  S. 

T  is  the  stage  transformation  function  that  expresses 
the  components  of  the  output  state  Y  as  a  function 
of  the  input  state  and  decisions,  Y=T(X,S). 

An  **'  after  an  Xi  or  Si  variable  indicates  the  optimal 
value  for  that  stage,  [Ref.  5  p22-23.  ]. 

The  one-stage  initial-state  optimization  problem  is  to  find 
the  best  stage  return  as  a  function  of  the  input  state  X. 
The  optimal  return  from  the  stage  will  be  denoted  as  F  (X) 
and  the  optimal  decision  as  S*  (X)  .  Thus, 

F(X)  =  min  R(X,S(X))=B(X,S*  (X)),  (2.6) 

S  (X) 

where  the  decision  S*{X)  provides  the  best  return  for  that 
particular  input  value  X.  If  X  were  to  change  the  values  of 
S*(X)  and  the  return  function  should  also  change. 

Now  consider  a  multistage  problem  as  shown  in  figure 
2.2.  This  is  the  standard  recursive  DP  structure  with  back¬ 
ward  numbering  where  stage  N  is  on  the  left  and  the  stage  1 
on  the  right.  The  optimal  return  for  stages  i,i-1,...,1  is 
represented  by  Fi(Xi).  The  solution  technigue  solves 
succeeding  stages  frcm  right  to  left  providing  the  optimal 
return  FN  (XN)  up  to  the  current  stage  under  consideration. 
In  the  single  stage  problem  Y  was  the  output  state  variable, 
tut  now  Y  serves  as  the  input  to  the  next  stage  so 
Xi- 1  =  Y=Ti (Xi,Si)  .  Thus,  in  an  N  stage  problem,  the  output 
from  stage  i  becomes  the  input  to  stage  i- 1  for  i=1,...,N. 

In  figure  2.2 

Xi  is  the  state  variable 
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Figure  2.2  Example  of  an  N  Stage  Decision  Problem. 


Si  is  the  decision  variable 
Ri  is  the  stage  return  function 
Ti  is  the  stage  transformation  function 
*  after  an  Xi  or  Si  variable  indicates  the  optimal 
value  for  that  stage# 

where  i  =  1,...,N  is  the  stage  index.  For  each  stage  a  deci¬ 
sion  Si*  (Xi)  will  be  made  that  yields  the  optimal  return 
Fi(Xi)  as  a  function  of  the  present  stage  return 
R  (Xi  #  Si*  (Xi) )  and  the  optimal  return  Fi-l(Xi-l)  for  the 
stages  1  through  i-1.  The  recursive  eguation  for  a  mini¬ 
mization  procedure  would  be# 

Fi  (Xi)  =min  (Ri  (Xi#  Si)  +Fi-  1  (Xi- 1) )  i  =  1,...,N#  (2.7) 

Si 

wher  e 

Xi- 1=7i  (Xi,S i)  =Xi-CiSi  i  =  1#...,N.  (2.8) 

The  stage  transformation  is  simply  the  budget  minus  the 
guantity  purchased  multiplied  by  its  cost.  This  recursion 
assumes  the  overall  return  function  is  the  weighted  sum  of 


the  individual  return  functions  for  all  replacement  parts. 
The  equation  for  this  relationship  is, 

MSRT=£Ei (Xi,Si)  i=1,...,N.  (2.9) 

The  problem  is  to  make  a  series  of  decisions  SN...S1  that 
will  optimice  the  overall  return  function  while  remaining 
within*  the  budget  constraints. 

For  the  previously  discussed  provisioning  problem,  using 
the  different  replacement  parts  as  the  stages,  the  state 
variable  Xi  is  the  resource,  i.e.,  the  dollars  remaining  at 
each  stage.  The  decision  variable  Si  is  the  number  of  parts 
purchased  in  stage  i.  The  stage  return  function  is  the  MSRT 
which  must  be  computed  from  four  operating  characteristics 
for  each  part,  (essentiality,  procurement  lead  time,  demand 
rate,  performance  as  a  function  of  quantity  procured).  The 
following  equation  represents  the  return  function  at  each 
stage. 

Hi  (Xi,Si)  =Ei*Li*Ti*Zi  (Si)  /_Ei*Li*Ti  (2.  10) 

1  =  1,. .  •  ,  N  . 

The  Standard  DP  approach  computes  the  function  Fi(Xi) 
from  equation  2.7  recursively  for  i  =  1,...,N  and  for  all 
possible  values  of  Xi.  At  each  stage  i  this  is  done  by 
examining  all  feasible  decisions  for  each  Xi,  computing  the 
return  for  each  decision,  and  then  adding  the  corresponding 
optimal  return  Fi(Xi-l)  from  the  remaining  stages  to  arrive 
at  the  optimal  return  for  the  present  Xi.  Then  for  all  Xi 
values  in  the  stage,  the  best  total  return  is  chosen  and 
assigned  to  Fi  (Xi)  .  At  the  Nth  stage  the  process  is 
complete  and  FN (XN)  represents  the  optimal  return  for  the 
problem.  This  process  is  efficient  for  small  problems, 
however  it  becomes  computationally  more  difficult  as  the 
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number  of  stages  and  the  number  of  possible  state  variable 
values  increase  in  each  stage. 

For  an  example  of  the  computational  burden,  consider  a 
very  simple  problem  with  a  budget  3,  a  constant  increment  K 
for  the  state  variable  Xi  within  each  stage,  and  a  max  of  D 
decisions  for  each  state  variable  value.  Ir.  this  case  the 
number  of  evaluations  for  one  stage  is  equal  to  (3/K)  *D. 
Assuming  the  transformation  function  and  the  recursive  equa¬ 
tion  for  the  optimal  return  equation  require  only  one  addi¬ 
tion  each,  then  there  are  2*(B/K)*D  additions  per  stage. 
Each  state  variable  value  will  require  D- 1  comparisons, 
hence,  there  are  (B/K)*(D-1)  comparisons  per  stage.  The 
total  number  of  operations  (additions  plus  comparisons) 
required  for  all  stages  may  be  represented  by. 


N  (  (2*B*D)/X)  *(B*  (D-1)/K)  , 


(2.  11) 


or  simplified  this  becomes 


(N  *B/K)  *  (3D-  1 )  . 


(2.  12) 


This  example  is  actually  a  gross  over-simplification  cf  the 
computational  problem.  It  is  possible  the  return  function 
could  be  very  complex  requiring  many  operations  for  each 
decision  at  each  Xi  value  through  the  range  of  the  budget. 
However,  for  a  numerical  comparison,  consider  a  25-stage 
problem  with  a  budget  of  74000,  a  state  variable  increment 
of  1000  and  a  maximum  of  100  decisions.  This  leads  to 
553150  operations  as  shown  in  equation  2.13, 


(  (2  5*74000) /I  0  00)  *  (  (3*  100  )  -  1 )  =5  53 , 1  50 . 


(2.  13) 


This  is  the  required  number  of  additions  and  comparisons  for 
a  25  stage  problem.  For  today's  high  speed  computer  an 
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addition  operation  at  machine  language  level  takes  about 
twenty-four  machine  cycles  while  a  comparison  takes 
eighteen.  Using  twenty  cycles  as  an  average  and  a  machine 
speed  of  two  million  cycles  per  second,  the  above  number  of 
comparisons  and  additions  would  take  approximately  nine 
seconds.  This  time  estimate  is  for  an  incredibly  simple 
problem  that  does  not  take  into  account  any  of  the  adminis¬ 
trative  overhead,  objective  or  constraint  computations  or 
increased  calculations  required  by  complex  functions.  3y 
setting  K=10,  which  is  equivalent  to  reducing  the  stage 
increment  from  1000  to  10,  equation  2.14  shows  that  the 
number  of  operations  will  increase  to 


{  (25*7  400  0) /I  0)  *  (  (3*100)  -  1)  =  55,  31  5,  00  0. 


(2.  14) 


Hence,  it  is  easy  to  see  the  limitation  of  standard  dynamic 
programming.  While  it  gives  integer  solutions,  the  computa¬ 
tional  burden  is  considerable  even  in  simple  problems. 


C.  A  DYNAMIC  PROGRAMMING  7 ARIATION 

The  computational  workload  of  the  full  dynamic  program¬ 
ming  procedure  makes  it  undesirable  for  large  problems.  A 
technique  that  will  be  referred  to  as  ' FUNNELING*  because  of 
its  control  of  the  resource  boundaries  between  the  stages, 
was  developed  and  applied  to  a  standard  recursive  D? 
computer  program.  The  result  is  a  shortened  procedure  that 
provides  a  locally  optimal  integer  solution  to  problems  that 
can  be  formulated  in  the  DP  construct. 

For  any  optimal  solution  to  an  N-stage  DP  problem  there 
is  a  series  of  decisions  Si*  for  i=1,...,N.  Corresponding 
to  this  solution  there  is  a  sequence  of  Xi*  values  for 
i=1,...,N.  These  optimal  state  variable  values  Xi*,  repre¬ 
sent  the  best  quantities  of  resource  available  for 


consumption  at  that  stage  and  are  not  known  in  advance.  If 
it  were  possible  to  find  an  approximate  value  for  Xi*  in 
each  stage,  the  optimal  solution  then  could  be  localized 
with  an  iterative  process.  That  is,  rather  than  compute 
F  (X)  for  all  X  values  there  could  be,  a  'window'  of 
interest,  centered  around  the  approximate  Xi*  value  in  each 
stage.  The  computations  would  be  performed  on  the  limited 
range  of  Xi  values  in  that  window  thereby  reducing  the 
required  workload  significantly.  This  concept  of  control¬ 
ling  computational  boundaries  around  an  approximate  optimal 
Xi*  value  is  the  foundation  for  the  D?  variation. 

The  problem  initially  is  how  to  localize  the  regions  of 
the  best  Xi.  This  is  solved  by  the  first  run  of  the  funnel 
program  using'  a  large  state  variable  increment.  The 
resulting  Xi  values  are  then  used  as  an  approximation  to  the 
initial  Xi*  values  i=1,...,N.  a  subroutine  of  the  funnel 
program  then  constructs  a  window  around  the  area  of  interest 
for  further  refinement.  The  first  estimate  of  the  Xi* 
values  is  quite  rough  and  the  decisions,  especially  the  Si 
which  have  costs  less  than  the  initializing  increment,  may 
cnange  significantly  as  the  Xi  are  examined  with  a  lesser 
state  variable  increment. 

The  second  program  run  is  the  first  refinement  of  the 
initial  solution.  Using  the  large  increments  the  program 
will  have  chosen  the  test  return  available  though  it  may  not 
be  close  to  the  optimal  value.  The  increased  sensitivity  of 
the  reduced  state  variable  increment  will  allow  considera¬ 
tion  of  more  Xi  values  in  the  window  permitting  the  opti¬ 
mization  process  to  move  closer  to  the  best  return.  It  is 
conceivable  that  any  Xi  may  change  as  much  as  the  window  of 
its  stage  will  allow.  This  much  variation  in  the  decisions 
indicates  that  the  best  solution  lies  outside  the  window  and 
the  boundaries  must  be  changed  accordingly.  This  situation 
may  be  present  in  several  stages  for  any  program  iteration. 
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Due  to  the  characteristics  of  the  DP  procedure, 
the  state  variable  consumed  at  stage  N  wiil  cause  the 
quantity  of  resource  available  at  stage  >!- 1  and  sur-ceeming 
stages  to  change.  The  concept  is  tnat  repetitive  runs  of 
the  funneling  program  will  allow  the  computat ionai  bounds 
for  each  stage  to  localize  and  center  upon  their  respective 
Xi*.  When  each  of  the  stages  has  chosen  its  Xi*,  the  Si* 
for  i=1,...,N  will  represent  the  best  solution  for  a  given 
objective  and  budget.  This  solution  will  be  a  local  and 
possibly  a  global  optimum.  Further  discussion  of  the 
program  specifics  may  be  found  in  Chap  3. 

Additional  computations  may  be  saved  in  the  funneling 
program  by  ranking  the  stages,  largest  to  smallest,  by  the 
cost  of  their  respective  parts.  Purchasing  the  most  expen¬ 
sive  items  first  will  reduce  the  magnitude  of  the  budget  and 
simplify  the  computations  in  the  remaining  stages.  The 
amount  of  computation  depends  on  three  things,  1)  the  width 
of  the  window  for  each  stage,  2)  the  number  of  state  vari¬ 
able  increments  within  the  window  and  3)  the  number  of  deci¬ 
sions  available  for  each  state  variable  increment. 

Consider  again  the  computational  example  of  section  3, 
where  the  number  of  operations  is  estimated  for  the  dynamic 
programming  solution  procedure.  To  use  that  analysis  for 
the  funnel  procedure,  B  in  equation  2.12  is  replaced  by  the 
width  of  the  window  W.  This  width  will  vary  between  stages 
but  this  example  will  assume  Wi  is  constant  and  equal  to 
10000.  Setting  the  remaining  factors  the  same,  K=10,  S=100, 
and  N=25,  the  number  of  operations  required  from  equation 
2.  1  2  is, 

{  (25*10  000) /I  0)  *  (  (3*  100)  -1)  =7,475,000.  (2.  15) 

This  is  roughly  13.57c  of  the  full  DP  requirement  computed  in 
equation  2.11,  a  marked  improvement  in  efficiency.  There 
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III.  THE  FUNNEL  PROGRAM 


This  chapter  discusses  the  operation  of  the  standard  DP 
program,  the  modifications  required  to  employ  the  funneling 
concept  and  some  of  the  implementation  problems  encountered 
in  the  process.  Specific  subroutine  operations  are 
explained  in  the  second  section. 

A.  MODIFICATION  FOR  THE  FUNNEL  PROCEDURE 

The  CP  variation,  called  the  Funnel  program,  is  a  modi¬ 
fied  version  of  a  standard  DP  computer  program  named  DPS. 
The  modified  program  will  be  referred  to  as  DP6  for  the  sake 
of  clarity  in  the  following  discussion.  The  DP 5  program  as 
described  in  [Ref.  3]  has  four  subroutines,  5T3P.E T,  STORE, 
TRANFM  and  DLIMIT.  The  following  list  describes  what  func¬ 
tion  each  subroutine  performs  in  the  DP  process. 

1)  STORE  -  enters  the  constants  to  be  used  in  the  other 

subroutines  and  the  main  program. 

2)  DLIMIT  -  defines  the  range  of  the  decision  values  Si 

for  any  value  of  Xi  in  stage  i. 

3)  STGRET  -  defines  and  computes  the  return  function  for 

any  values  of  Xi  and  Si  in  stage  i,  Ri(Xi,Si). 

4)  TRANFM  -  defines  and  computes  the  transformation 

function  for  any  values  of  Xi  and  Si  in  stage 
i,  Ti  { Xi,  Si)  . 

Providing  DP5  with  the  specific  input  data  as  required  in 
[Ref.  3],  the  program  will  compute  the  optimal  return  by  the 
standard  recursive  DP  procedure.  The  following  figure  shows 
the  flow  diagram  for  the  general  solution  technique. 


Figure  3. 
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Compute  recursively 
Fi (X)=min(Ri(X,S)  +  Fi-l)  ! 

Si  I 

for  1*1,. ...N 

’  I 

using  DLIXIT  ; 

TRANS FM  j 

STCRET  ! 


1  Flow  Diagram  of  the  DP5  Procedure 
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The  input  data  are  the  parameters  for  the  particular  problem 
being  solved.  When  DP5  is  run  it  uses  the  STOPS  subroutine 
to  place  the  parameters  in  memory,  then  calls  the  other 
subroutines,  DLIMIT,  TRANFM,  and  STSRET  to  perform  the  DP 
procedure. 

To  use  the  concept  of  funneling  described  in  Chap  2, 
additional  subroutines  must  be  added  to  DP5  to  control  the 
state  variable  computational  boundaries  in  each  of  the 
stages.  There  are  two  possible  starting  conditions  for  the 
funnel  procedure,  1)  some  input  information  is  available  for 
the  decisions  at  each  stage,  or  2)  no  prior  knowledge  of  the 
decisions  is  available.  To  respond  to  these  conditions  two 
subroutines  ENTER  and  ADJUST  have  been  created  to  correctly 
manipulate  the  inputs.  The  following  list  describes  what 
function  each  of  the  subroutines  performs  in  the  Funnel 
program. 

1)  ENTER  -  reads  a  decision  set  from  a  file  then 

creates  the  input  data  file  in  the.  correct 
format  for  the  Funnel  program, 

2)  ADJUST  -  manipulates  the  computational  boundaries  of 

an  input  data  file  to  reflect  the  changes 
required  at  each  program  iteration. 

The  general  procedure  for  DP6  is  shown  in  figure  3.2, 

The  input  parameters  are  read  into  memory  with  the  STORE 
function.  Then,  depending  on  the  starting  condition,  the 
ENTER  function  may  used  to  create  an  input  data  file  for  the 
Funnel  program  from  a  previously  available  solution.  This 
solution  may  be  obtained  from,  1)  some  other  solution 
method,  2)  a  rule-of-thumb  estimate  of  the  solution,  or  from 

3)  past  information  on  similar  problems.  In  any  of  these 
cases  the  estimated  or  best  guess  decision  set  is  trans¬ 
formed  to  the  correct  input  format  for  the  Funnel  program. 
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When  there  is  no  previous  solution,  the  beginning  data  rile 
may  be  constructed  in  the  correct  format  as  defined  by  the 
user  instructions  located  in  the  comment  section  of  DP 6,  see 
appendix.  This  input  file,  as  described,  may  be  used 
directly  with  the  program.  The  next  step  in  the  figure  is 
the  first  iteration  of  the  Funnel  process  utilizing  the 
DLI M IT,  TRANFM,  and  STGRET  subroutines.  When  the  calcula¬ 
tions  are  complete  the  results  must  be  examined  to  see  if 
the  process  has  found  the  optimal  solution,  i.e.,  ’is  it 
DONE’.  The  requirements  for  the  'DONE'  condition  are  1}  an 
all  integer  solution,  2)  no  window  boundaries  violated  and 
3)  the  budget  remaining  must  be  less  than  the  smallest  cost 
for  items  to  be  purchased.  If  all  of  these  conditions  are 
not  met  the  ADJUST  subroutine  is  activated  and  the  input 
data  file  is  manipulated  to  correct  the  problem  areas.  With 
each  program  iteration  this  corrective  procedura  is  repeated 
until  all  the  conditions  are  met  and  an  optimal  solution  is 
produced.  An  example  of  the  Funnel  program  output  is  shown 
is  figure  3.3. 


THE  OBJECTIVE  IS  MSRT 
CONST HAINT  ON  COST 

THB  DESCRIPTION  WHICH  APPEARS  BEIC1  FOR  STAGE  1 
APPLIES  TG  STAGE  1  THRU  STAGE  3  INCLUSIVE 
THE  PROBLEM  IS  TO  MINIMIZE  A  3  STAGE  PROCESSES 
IS  TO  BE  CHOOSEN  OPTIMALLY  BETWEEN  XN  =  2- 05000E*02 
AND  XN=2.Q5QQD*02 

OPTIMAL  XN=2.05000D*02  OPTIMAL  BETDRN=5-30246D*00 
N  XN  SN  XN-1 

3  2. 0500D*02  1.00D*01  5.50000*01 

2  5.50000*01  2. OOD*CO  3.50C0E*01 

1  3- 5000D*0 1  7.00D*00  0.0 


Figure  3.3  Example  of  DP6  Output 


The  objective  and  constraint  are  stated  in  the  first  two 
lines.  The  description  that  follows  gives  the  number  of 
stages,  the  budget  quantity,  the  optimal  return  and  detailed 
stage  results.  Specifically  delineated  for  each  of  the 
stages  are  the  input  budget  quantities,  the  decisions  and 
the  output  budget  quantities.  For  this  particular  example 
the  budget  remaining  at  the  end  of  the  problem  is  zero. 

Two  ancillary  functions  were  added,  one  to  the  main 
program  and  one  to  a  subroutine  to  provide  information  to 
the  user.  The  first  is  a  timing  function  that  provides  the 
program  run  time.  The  second  is  a  counting  function  that 
presents  to  the  screen  and  stores  in  a  file  the  numbers  of 
the  stages  in  which  the  optimal  return  values  are  too  close 
to  the  computational  boundaries.  Both  functions  have 
provided  valuable  insights  that  will  be  discussed  in  Chapter 
5. 

B.  THE  SUBROUTINE  ADJUST 

The  subroutine  ADJUST  is  the  heart  of  the  tunneling  tech¬ 
nique.  To  gain  the  computation  reduction  desired,  this 
subroutine  manipulates  the  windows  around  the  ’current’ 
state  variable  values,  Xi’  ,  in  each  stage  for  each  program 
iteration.  Initially  the  upper  and  lower  limits  HXi’  and 
Lxi'  on  the  window  at  stage  i  are  defined  by, 

HXi*  =  Xi'  +K*Ci+  1 ,  (3.1) 


and 


L  Xi '  =  Xi'-K*Ci+1,  (3.2) 

where  K  is  a  constant.  These  bounds  are  illustrated  in  fig 
3.4.  Further  adjustment  of  these  limits  may  be  necessary  to 
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assure  that  one  of  the  Xi  values  considered  in  the  computa¬ 
tion  of  F  (Xi)  is  Xi’.  At  stage  i  in  the  computation  of 
F(Xi)  the  optimization  considers  first  Xi=LXi/  then  Xi  is 
incremented  by  Ci  to  give  the  next  value,  etc.  Thus  if  IXi 


is  not  properly  set,  the  value  Xi=Xi’  will  never  be 
considered. 


Figure  3.4  Example  of  a  iindow  on  a  Stage  i 


This  adjustment  of  LXi'  is  made  by  decrementing  from  the 
temporary  stage  value,  Xi'  ,  by  the  cost  Ci  for  that  stage. 
When  this  decrementing  process  first  yields  a  value  less 
than  or  equal  to  LXi'  that  value  is  set  as  the  nt*  lower 
limit  on  the  window  and  is  called  LXi.  A  similar  adjustment 
is  made  to  the  upper  limit  by  incrementing  from  Xi'  until 
the  value  HXi'  is  equalled  or  exceeded.  The  new  upper  limit 
is  called  HXi.  This  process  is  repeated  for  ail  the  stages 
until  the  funnel  is  formed.  In  figure  3.5,  tne  funnel  may 
be  seen  for  N  stages  of  a  provisioning  problem.  The  budget 
quantity  3  is  represented  by  the  vertical  lines  for  all 
stages.  The  upper  and  lower  bounds  are  represented  by  HXi, 
LXi  for  i=1,...,N  where  i  is  the  stage  index.  In  each  stage 
the  bounds  are  at  least  a  distance  equal  to  the  product  of 
the  constant  K  and  Ci+1  from  the  Xi*  value.  Stages  where 
the  bounds  are  further  from  the  center  indicate  the  addi¬ 
tional  adjustment  for  compatability  between  the  successive 
Xi*  values. 


Lgure  3.5  A  Funnel  Diagram  for  N  Stages  of  a  Problem 


> 


The  reason  for  multiplying  the  succeeding  stage  costs  by 
the  constant  K  is  to  allow  the  DP  procedure  some  computa¬ 
tional  freedom  in  each  of  the  stages  while  searching  for  the 
Xi  and  3i*(Xi)  values.  Daring  the  computational  process, 
should  a  stage  choose  an  Xi'  tnat  is  within  Ci+1  of  the 
computational  bound,  the  boundary  is  said  to  be  violated. 
This  is  because  the  DP  procedure  may  actually  prefer  a  value 
outside  the  bound  but  is  unable  to  reach  it.  Figures  3.6 
through  figure  3.9  show  an  example  of  the  stage  boundary 
manipulation  as  the  iterative  process  of  the  DPb  program 
searches  for  the  optimal  solution. 


—  200  =hX  i 


-150  =  X  i  ' 


-100  =LXi 


i-l 

-j-170  =HXi-  1 

—  L2Q.-Xt-1  ' 

,170  =  LX  i  -  1 


i-2 


120  =  H  X  i  -  2 


70  =  X  i  -  2  ’ 


20  =  L X  i  -  2 


Figure  3.6  The  Initialization  Run. 


The  initializing  run  depends  on  the  starting  condition  for 
the  problem.  If  the  ENTER  function  is  used,  the  state  vari¬ 
able  increments  will  represent  the  costs  of  the  respective 
stages.  Figure  3.6  shows  a  starting  condition  of  no  prior 
knowledge;  hence,  the  state  variable  increment  is  a 
constant,  in  this  example  10,  for  all  stages  and  of  egual 
number  on  either  side  of  the  current  Xi'  estimate. 

A  second  run  of  the  program  provides  an  opportunity  for  the 
DP  process  to  evaluate  the  estimated  optimal  returns  with  a 


I 
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i  2  00  =  HX i 


i-  1 

_ 1 70  =  HXi -  1 


U0=Xi  ' 

■  100=LXi 
Ci=  10 


i-2 

120=HXi-2 


=  X i-  1  ' 
d_7Q=LXi = 1 _ 


Ci-  1  =  5 


S0  =  Xi-2  ' 


±  20=LXi-2 
C i -  2  =  4 


Figure  3.7  The  Second  Run. 


smaller  state  variable  increment  or  a  revised  boundary  value 
and  to  adjust  the  window  around  a  new  Xi '  accordingly. 
Figure  3.7  shows  the  funnel  after  the  second  run  tut  before 
the  windows  are  adjusted.  The  program  has  selected  Xi'=l40, 
Xi-1' =70,  and  Xi-2'=50.  These  are  the  current  estimates  of 
the  optimal  Xi  •  values  at  each  stage.  Notice  that  stage  i- 1 
has  a  boundary  that  is  violated,  i.e.,  the  Xi-1 '  for  the 
stage  is  within  Ci=10  of  the  boundary.  In  this  case  the 
optimization  function  may  prefer  an  Xi-1  outside  the  window 
but  the  boundary  prevents  its  selection.  Also  in  the  second 
run  the  state  variable  increments  have  been  changed  to  the 
cost  of  the  item  associated  with  the  stage.  In  most  stages 
these  values  are  smaller  then  the  initial  increment  making 
the  return  functions  more  sensitive.  This  increase  in 
sensitivity  should  improve  the  optimization  selection. 

The  violation  of  the  lower  boundary  in  stage  i- 1  causes 
the  iteration  to  fail  the  conditions  of  the  DONE  checkpoint 
of  figure  3.2.  Conseguentl y ,  the  ADJUST  subroutine  is  acti¬ 
vated  again  to  manipulate  the  boundaries  producing  the 
results  shown  in  figure  3.8. 
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r  1  9  0  =  H  X  i 
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1 20=HXi- 1 
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L  1  4  0  =  X i  ’ 
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• 74=HXi-2 
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f  5  0  =  X  i  -  2  ' 

-  9  0=  L X  i 

f 

I 
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l 

20=LXi-l 

I 

1  1  8  =  L  X  i  -  2 

c  i  = : 

I  0 

Ci-l 

=  5 

Ci- 2  =  4 

Figure  3.8  The  Window  Adjustments. 

The  points  in  this  figure  represent  the  input  file  for  the 
third  run  of  the  program.  In  stage  i- 1  the  window  bound¬ 
aries  are  centered  cn  the  previous  iteration  LXi-1  value. 
Notice  the  window  bounds  for  stage  i-2  have  automatically 
adjusted  downward  to  compensate  for  the  increased  resource 
consumption  and  avoid  cascading  boundary  violations. 


r  1  9  0  =  H  X  i 

t 

I  1  4  0  =  X  i  ' 


9  0  =  I.  X  i 


|120=HXi-  1 

T 

|  65  =  Xi-  1  ' 

] 

1  2  0=  LX i -  1 


-  74=HXi-2 

i  4  6  =  X  i  -  2  ’ 
►  - 

1  - 
^ 18=LXi-2 


Ci=  10 


C i-  1=  5 


C  i  -  2  =  4 


Figure  3.9  The  Third  Run. 


Figure  3.9  shows  the  funnel  after  the  third  run.  The 
program  has  selected  Xi'=140,  Xi-1’=65  and  Xi-2'=46  for  the 
current  estimates  of  the  Xi *  values  at  each  stage.  Notice 
that  the  Xi-1  '  value  is  just  outside  tne  previous  LXi-1, 
which  was  70.  In  this  example  the  Funnel  program  manipu¬ 
lated  the  boundaries  of  the  windows  to  allow  the  DP  process 
to  choose  different  Xi 1  values  for  the  three  stages,  i,  i-1, 
i-2.  Subsequent  runs  of  the  program  will  determine  the  Xi* 
for  each  stage.  The  process  is  'DONE*  when  the  remaining 
budget  <=  minimum  Ci,  an  optimal  integer  solution  is 
provided,  and  there  are  no  boundaries  violated. 

C.  I IMPLEMENTATION  PBOBLEHS 

From  the  characteristics  of  the  funneling  process,  there 
are  problems  that  arise  when  transitioning  between  stages. 
When  computing  the  maximum  and  minimum  value  of  the  decision 
to  be  made  for  each  Xi  in  a  stage  i,  the  upper  and  lower 
bounds  and  item  cost  cf  stage  i-1  must  be  considered. 

In  figure  3.10,  the  left  side  is  the  window  for  stage  i,  the 
right  side  is  the  window  for  stage  i-1.  For  the  stages  i, 
i-1,  the  highest  values  of  the  state  variables  are  repre¬ 
sented  by  HXi  and  HXi-1  while  the  lowest  values  are  repre¬ 
sented  by  LXi,  LXi-1.  In  the  program  DP6,  the  two  variables 
DIO W  and  DHIGH  represent  the  minimum  and  maximum  decisions 
available  for  the  state  variable  in  each  stage.  It  is 
essential  that  DHIGH  for  any  Xi  value  in  the  window  of  the 
ith  stage  be  of  such  a  magnitude  that  it  will  not  force  the 
Xi-1  value  below  the  lower  bound  LXi-1.  Additionally,  DLOW 
must  be  of  such  magnitude  that  for  any  Xi  in  the  window  of 
the  ith  stage,  it  will  not  force  the  Xi-1  value  to  violate 
the  upper  bound  HXi-1.  If  either  of  taese  decision  quanti¬ 
ties  is  incorrect  at  stage  i  the  resource  quantity  will  be 
outside  the  window  at  stage  i-1.  The  derivation  of  the 
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Figure  3.10  Boundary  Manipulation. 


formulas  for  these  values  follows.  The  transformation  func¬ 
tion  as  represented  by  equation  2.8  of  Chap  2  is. 


Xi-1=Xi-CiSi  i  =  1,...,N. 


(3.3) 


If  this  transformation  is  to  take  place  from  a  window  in 
stage  i  to  a  window  in  stage  i-1  the  following  conditions 
must  be  true. 


LXi-1  <=  Xi- Ci* S  HIGH, 


(3.4) 


where  SHIGH  is  the  maximum  number  of  items  to  be  purchased 


X i-Ci *SLCW  <=  HXi-1, 


(3.5) 


where  SLOW  is  the  minimum  number  of  parts  to  be  purchased. 
By  algebraic  manipulation  a  relationship  for  SLOW  and  SHIGH 
may  be  established. 


SHIGH  <=  (Xi-LXi-1)/Ci, 


(3.6) 


and 

SLOW  >=  (Xi-HXi-  1)  /Ci.  (3.7) 

The  expressions  SHIGH  AND  SLOW  are  equivalent  to  the  DHIGH, 
DLDW  required  in  DP6  and  are  coded  in  the  subroutine  D LI M IT 
to  control  the  boundary  calculations.  It  is  conceivable 
that  at  the  extreme  stare  variable  values  of  Xi  the  computed 
boundaries  could  become  negative  or  larger  than  the  assigned 
budget.  Thus,  the  decisons  are  further  limited  by  DHIGH  = 
min  (SHIGH, 200)  and  DICW=  nax(0,  SLOW),  where  200  is  an  arbi¬ 
trarily  chosen  value. 

The  initialization  for  the  problem  is  either  from  input 
decisions  or  else  a  large  state  variable  increment  is  used. 
After  the  subroutine  ADJUST  manipulates  tne  windows  the 
first  time,  the  state  variable  increments  are  decreased  or 
the  boundary  conditions  are  changed  to  allow  examination  of 
the  new  possible  Xi  values  for  the  best  return.  The  subrou¬ 
tine  DLINIT  automatically  adjusts  the  DHIGH  and  DLOW  values 
to  compensate  for  the  changes  in  the  state  variable  incre¬ 
ment.  With  this  closer  scrutiny,  the  optimal  return  value 
may  change.  The  ADJUST  subroutine  moves  the  boundaries  up 
or  down  accordingly  when  an  optimal  value  is  within  one 
increment  of  the  bounds.  Eecall  that  if  the  optimal  return 
value  is  within  one  increment  of  window  bound,  it  may  be 
that  the  D?  process  would  prefer  a  value  outside  the  window. 
Hence,  the  boundary  is  recorded  as  violated,  and  the  window 
is  moved  for  the  next  iteration.  This  process  is  repeated 
until  all  the  optimal  return  values  are  centered  within  the 
windows  for  all  stages. 

To  save  computations  and  avoid  computing  Xi  values  that 
are  not  compatible  between  stages,  the  state  variable 


increment  for  each  stage  is  the  cost  of  the  item  to  be 
purchased  in  the  previous  stage.  A  problem  arises  in  that 
the  original  DP5  program  may  provide  solutions  that  are  not 
integer.  This  is  possible  because  in  the  stage  being 
considered,  the  best  return  may  reguire  an  Xi •  value  from 
succeeding  stages  that  has  not  been  computed.  DP5  handles 
this  by  interpolating  linearly  between  the  state  variable 
values  it  does  have  to  arrive  at  an  estimate  for  the  one 
required.  In  DP6  tc  avoid  this  approximation  and  obtain 
integer  solutions,  the  subroutine  ADJUST,  (as  discussed  in 
section  A),  manipulates  the  lower  and  upper  window  values  to 
ensure  the  previous  stage  optimal  return  values  are  avail¬ 
able  for  the  Xi  being  considered.  This  adjustment  process 
is  also  repetitive  and  is  done  an  conjunction  with  the 
boundary  control. 

Another  problem  arises  from  the  effort  to  save  computer 
computation  time.  If  there  is  a  wide  range  of  costs,  the 
use  of  the  item  cost  as  the  state  variable  increment  within 
a  stage  becomes  troublesome.  Initially  the  ADJUST  subrou¬ 
tine  creates  the  state  variable  windows  for  the  individual 
stages  using  a  constant  multiplier  that  may  be  set  by  the 
user.  To  remedy  the  problem  of  the  extreme  range  of  the 
costs  of  items  for  the  provisioning  problem  being  consid¬ 
ered,  ($10000.00  to  .50),  the  multiplier  was  changed  to  a 
variable  that  depends  directly  on  the  item  cost.  This  was 
done  to  decrease  the  program  repetitions  required  to  reach 
an  optimal  solution.  Under  the  constant  multiplier  method, 
a  change  in  the  purchase  quantity  of  a  larger  cost  item, 
(i.e.$1500),  could  require  many  program  repetitions  to 
re-allocate  the  resource  available  for  the  lower  cost  items, 
(i. e. $10,5, 3.75,  etc).  The  smaller  window  range,  due  tc  the 
constant  multiplier,  hampers  the  re-allocation  effort.  The 
variable  multiplier  allows  greater  compensation  for  lower 
cost  items  to  facilitate  a  more  rapid  change  if  required. 
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17.  ANALISI S 


h 


I 

i. 
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The  primary  topic  of  this  chapter  will  be  the  comparison 
of  the  Marginal  Analysis  and  Tunneling  solution  techniques 
for  the  previously  stated  problem.  The  discussion  will 
cover  the  results  of  both  programs  and  the  solution  differ¬ 
ence.  The  standard  dynamic  program  procedure  will  not  be 
addressed  because  of  the  computational  advantage  of  the  DP 
variant  and  the  similiarity  of  the  two  techniques. 

Tor  simplicity  and  ease  of  comparison  a  small  three 
stage  problem  will  be  defined  and  solved  using  both  solution 
procedures.  As  in  the  original  provisioning  problem  the 
objective  is  to  minimize  the  MSRT  with  a  constraint  on  cost. 
Three  different  budget  quantities  will  be  examined  to  demon¬ 
strate  the  differences  in  the  solution  techniques.  The 
following  table  is  a  list  of  the  input  parameters  which  are 
cost  Ci,  demand  rate  Li,  essentiality  Ei,  and  procurement 


lead  time  Ti. 

For  this 

particular 

example  Ti  is  one  year. 

TABLE  II 

Sample  Problem  Input 

Parameters 

I 

DEMAND 

COST 

T  ESSEN 

1 

5.000 

5.00 

1.00  3.0 

2 

1.000 

10.00 

1.00  2.0 

3 

10.000 

15.00 

1.00  1.0 
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For  the  first  budget  quantity  of  3195.00  the  Marginal 
Analysis  program  returns  the  following  output. 


TABLE  III 

Marginal  Analysis  Program  Return,  Budget=$ 195. 00 

*********************************************** 

*  METHOD:  PIXED  BUDGET  =  195.00  * 

*********************************************** 


********************************** 


STEP 

ITEM 

SI 

BATIO 

1 

1 

1 

0.480808496 

2 

1 

2 

0.365659833 

3 

1 

3 

0.260618091 

4 

1 

4 

0.  172421157 

5 

1 

5 

0.  105280340 

6 

2 

1 

0.0735758543 

7 

3 

1 

0.0600002967 

8 

1 

6 

0.0591956675 

9 

3 

2 

0.0533369593 

10 

3 

3 

0.046688754 1 

11 

3 

4 

0.0400909968 

12 

3 

5 

0.0336193480 

13 

1 

7 

0.0306577347 

14 

3 

6 

0.0273999237 

15 

3 

7 

0.0216008648 

16 

2 

2  • 

0.0207276568 

17 

3 

a 

0.0164023377 

18 

1 

8 

0.0146531016 

19 

3 

9 

0.0119544677 

****** 

BB 

190. 

185. 

180. 

175. 

170. 

160. 

145. 

140. 

125. 

110. 

95. 

80. 

75. 

60. 

45. 

35. 

20. 

15. 

0. 


******* 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

0 


*********************************************** 

MINIMIZE  HSET  EXCESS:  0.0 

I  SI  MS8T 

1  8  1.302 

2  2  10.396 

3  9  12.176 

OVEBALL  6.003 

REALIZED  PEBPORM ANCE  =  6.0033 

BOUNDS  (  6.0033,  6.0033) 


The  decision  results  are,  S 1=8,  S2=2,  S3=9  and  the  return  is 
Z=  6.0033,  with  no  budget  remaining.  Viewing  the  purchase 
sequence  in  table  III  it  may  be  seen  that  the  Marginal 
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Analysis  program  has  reached  a  solution  uninterrupted,  that 
is  without  implementing  the  heuristic  developed  in  [Ref.  1]. 
Consequently,  for  this  budget  quantity  tne  solution  is 
optimal.  Notice  another  indicator  of  the  optimal  solution 
in  the  last  line  of  table  III.  The  bounds  described  are  the 
possible  upper  and  lower  values  for  the  return.  When  these 
quantities  are  the  same  value,  the  solution  is  the  best 
available. 


i 


i 

The  same  budget  and  parameters  entered  in  the  Funnel 
program  produce  the  results  in  table  IV.  The  output  is  in 
the  same  format  discussed  in  Chapter  3.  The  initial  budget 
1  quantity,  XN  and  the  optimal  return  are  shown  in  the  line 

just  above  the  individual  stage  results.  The  variable  XI  is 
the  budget  quantity  at  the  beginning  of  the  stage.  Si  is  the 
decision  or  quantity  purchased,  and  Xi-1  is  the  budget  at 
the  end  cf  the  stage.  Recall  that  item  3  is  the  most  expen¬ 
sive  and  item  1  the  least  expensive.  The  optimal  solution 
foe  the  $195.00  budget  is  S1*  =  8,  S2*=2,  and  S3*  =  9  with  a  Z 

value  equal  to  6.0033  and  no  budget  remaining.  This  answer 


TABLE  If 

Funnel  Program  Results,  Budget=$ 195. 00 

THE  OBJECTIVE  IS  BSBT 
CONSTfi AIST  OM  COST 

THE  DESCRIPTION  WHICH  APPEABS  BE1CW  FOB  STAGE  1 
APPLIES  TO  STAGE  1  THRU  STAGE  3  INCLUSIVE 
THE  PROBLEM  IS  TO  MINIMIZE  A  3  STAGE  PROCESSXN 
IS  TO  BE  CHOOSER  OPTIMALLY  BETWEEN  X N  = 1 . 9 5 00 D *0 2 
AND  XN=  1. 9500D  +  02 

OPTIMAL  XN=1.95Q0QD*Q2  OPTIMAL  EETD8N=6.Q0332D*00 
N  XN  SN  XN—  1 

3  1.9500D*02  9. Q0D*  00  6.0000D*01 

2  6. 0000  D*  0 1  2.00D*00  4.0000D*01 

1  4. OOOOD+  0 1  8.00D  +  00  0-0 


5 


v  'it*  v 


matches  exactly  the  Marginal  Analysis  program  results  shown 
in  table  III,  confirming  tnat  the  Marginal  Analysis  solution 
is  optimal. 


TABLE  V 

Marginal  Analysis  Program  Return,  Budget=$200. 00 

- ***************  **  ****  *s  ***  ******  ** '* ******  ****** 

*  METHOD :  PIXED  EODGET  =  200.00  * 

*********************************************** 

*********************************************** 


STEP 

ITEM 

SI 

RATIO 

BR 

1 

1 

1 

0.480808496 

195.00 

2 

1 

2 

0.365659833 

190.00 

3 

1 

3 

0.260618091 

185.00 

4 

1 

4 

0.  172421157 

180.00 

5 

1 

5 

0.105280340 

175.00 

6 

2 

1 

0.0735758543 

165.00 

7 

3 

1 

0.0600002967 

150.00 

8 

1 

6 

0.0591956675 

145.00 

9 

3 

2 

0.0533369593 

130.00 

10 

3 

3 

0.0466887541 

115.00 

11 

3 

4 

0.0400909968 

100.00 

12 

3 

5 

0.0336193480 

85.00 

13 

1 

7 

0-0306577347 

80.00 

14 

3 

6 

0.0273999237 

65.00 

15 

3 

7 

0.0216008648 

50.  00 

16 

2 

2 

0. 02072765t>8 

40.00 

17 

3 

8 

0.0164023377 

25.00 

18 

1 

8 

0.0146531016 

20.00 

19 

3 

9 

0.0119544677 

5.00 

20 

3 

10 

0.0083406679 

-10.00 

20 

1 

9 

0.0064818673 

0.0 

*********************************************** 

MINIMIZE  MSHT  EXCESS:  0.0 


1 
1 

2 
3 


SI 

9 

2 

9 


MSB! 
0.514 
10.396 
12. 176 


OVERALL  5.565 


REALIZED  PERFORMANCE 
BOUNDS  (  6.0033, 


=  5.565 

4.3120) 


For  the  second  example  the  budget  quantity  of  $200.00 
was  used.  The  Marginal  Analysis  program  returns  the  output 


in  table  V.  The  results  are  S3  =  9,  52=2,  and  S1=9  with  a  Z 
value  of  5.5652  and  no  budget  left  over.  Examining  the 
detailed  output  of  the  Marginal  Analysis  program,  table  V, 
it  may  be  seen  that  the  purchase  of  the  twentieth  item 
drives  the  budget  negative.  The  program  implements  the 
heuristic,  re-evaluating  the  benefit  to  cost  ratios,  then 
chooses  the  number  one  item  for  the  next  purchase.  Notice 
that  in  the  bottom  line  'of  table  7  the  upper  and  lower 
bounds  are  not  equal.  This  is  a  good  indication  the  solu¬ 
tion  may  not  be  optimal  and  should  be  verified. 


TABLE  71 

Funnel  Prograa  Results,  Budget=$200. 00 


THE  OBJECTIVE  IS  MSHT 
CONSTRAINT  ON  COST 

TBE  DESCRIPTION  BHICH  APPEARS  EEICW  FOR  STAG!  1 
APPLIES  TC  STAGE  1  THRU  STAGE  3  INCLUSIVE 
THE  PROBLEM  IS  TO  MINIMIZE  A  3  STAGE  PBOCESSXN 
IS  TO  BE  CHOOSEN  OPTIMALLY  BETHEEN  XN=2. 00C 0D *02 
AND  XN=2. 000QD+Q2 

OPTIMAL  XN=2. 0000D+02  OPTIMAL  R£T0RN=5. 565 19D*00 
N  XN  SN  XN-1 

3  2. 0000D  +  02  9,Q0D#-00  6.500CD*01 

2  6. 500  0D+  0 1  2. 0  CD  +  CQ  4.5000D+01 

1  4.5000D*01  9. 00D+Q0  0-0 


The  same  budget  and  parameters  entered  in  the  Funnel 
program  produce  the  results  given  in  table  VI.  In  this  case 
the  answer  is  again  identical,  S1=9,  S2=2,  and  S3=9  with  the 
optimal  return  Z=  5.  56519.  Incidentally,  the  heuristic  in 
the  Marginal  Analysis  program  makes  the  correct  choice  at 
the  twentieth  step  when  it  selects  the  item  with  the  second 
highest  benefit  to  cost  ratio. 

For  a  third  case,  a  budget  of  $205.00  is  entered  in  both 
programs.  The  results  for  the  Marginal  Analysis  code  are 
shown  in  table  VII. 


TABLE  TII 


Marginal  Analysis  Program  Hesults,  Budget=$205. 00 

*********************************************** 

*  METHOD:  FIXED  EUDGET  =  205.00  * 

*********************************************** 

*********************************************** 


STEP 

ITEM 

SI 

BATIO 

BB 

1 

1 

1 

0-480808496 

200.00 

2 

1 

.  2 

0.365659833 

195.00 

3 

1 

3 

0.260618091 

190.00 

4 

1 

4 

0. 172421157 

185.00 

5 

1 

5 

0. 105280340 

180.00 

6 

2 

1 

0.0735753543 

170.00 

7 

3 

1 

0.0600002967 

155.00 

8 

1 

6 

0.0591956675 

150.00 

9 

3 

2 

0.0533369593 

135.00 

10 

3 

3 

0.046683754  1 

120.00 

11 

3 

4 

Q.04C09 09968 

105. 00 

12 

3 

5 

0.0336193480 

90.00 

13 

1 

7 

0.0306577347 

85.00 

14 

3 

6 

0.0273999237 

70.00 

15 

3 

7 

0.0216008648 

55.00 

16 

2 

2 

0.0207276568 

45.00 

17 

3 

8 

0.0164023377 

30.  00 

18 

1 

8 

0.0146531016 

25.00 

19 

3 

9 

0.0119544677 

10.00 

20 

3 

10 

0.0083406679 

-5.00 

20 

1 

9 

0.0064818673 

5.00 

21 

1 

10 

0.0026624922 

0.0 

********* ** *********** ************************* 

MINIMIZE  BSHT  EXCESS:  0.0 

I  SI  HSBT 

1  10  0.190 

2  2  10.396 

3  9  12.176 

OVEBALL  5.385 

HEALIZBD  PEBFGBM ANCE  =  5-3852 

BOUNDS  (  o.  0033 ,  <*.3120) 


The  solution  is  S1  =  10,  S2=2,  and  S3=9  with  a  return  value  of 
Z=  5.3852.  The  upper  and  lower  bounds  are  not  equal  indi¬ 
cating  the  heuristic  has  been  applied  and  the  solution 
should  be  verified.  Examining  the  output  from  the  Marginal 
Analysis  program,  table  VI  I,  it  nay  be  seen  the  budget  was 
again  driven  negative  at  the  twentieth  item.  The  heuristic 
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was 


activated  purchasing  two  of  the  number  one  items  to 
drive  the  budget  to  zero. 

Entering  the  same  input  and  budget  parameters  in  the 
Funnel  program  produces  the  results  shown  in  table  VIII. 


TABLE  VIII 

Funnel  Program  Results,  Budget=$205. 00 

THE  OBJECTIVE  IS  MSRT 
CONSTRAINT  ON  COST 

THB  DESCRIPTION  WHICH  APPEARS  BEICN  FOR  STAGE  1 
APPLIES  TC  STAGE  1  THRU  STAGE  3  INCLUSIVE 
THE  PROBLEM  IS  TO  MINIMIZE  A  3  STAGE  PBOCESSXN 
IS  TO  BE  CHOOSEN  OPTIMALLY  EET8EEN  XN=2- 05000D+02 
AND  XN=2.Q500D*02 

OPTIMAL  XN=2.05000D*02  OPTIMAL  BET0RN=5. 30246D*00 
N  XN  SN  XN-1 

3  2-05000*02  1 . 00D*  0 1  5-5000D*01 

2  5.5000D*01  2. Q0D*C0  3.50C0D*01 

1  3- 500  OD*  0 1  7-00D*C0  0.0 


The  answers  are  not  the  same.  The  Funnel  program  results 
indicate  the  best  part  selection  is  S1  =  7,  52=2,  and  S3=10 

for  an  optimal  return  of  Z=  5.30246.  This  selection 

provides  an  MS RT  that  is  roughly  .08  better  than  the 

marginal  analysis  method. 

Recall  from  section  A  of  Chapter  2  that  the  Marginal 
Analysis  program  has  a  function  called  KIT  that  will  deter¬ 
mine  the  return  from  a  given  decision  set.  Entering  the 
optimal  solution  from  the  Funnel  program  into  the  KIT  func¬ 
tion  produces  the  result  in  table  IX.  The  resulting  objec¬ 
tive  value  is  identical  to  the  value  computed  by  the  DP 

variant.  This  shows  that  the  computational  procedures  are 
the  same  for  both  programs. 

Now  consider  the  definition  of  a  much  larger  problem. 
It  is  of  the  same  type  as  the  previous  examples,  an 


TABLE  IX 

Kit  Function  Results,  8 ud get=$205. 00 


EVALUATION  OF  SPECIFIED  ALLOCATION. .. . 

I  SI  BSRT 

1  7  3.085 

2  2  10.396 

3  10  7.610 

*  ************ *2 v**Ai £****;•***♦****•*♦» **•*•« 


objective  of  minimizing  MSP.T  with  a  constraint  on  cost.  The 
input  parameters  cost  Ci,  demand  cate  Li,  essentiality  Ei, 
and  procurement  lead  time  Ti,  are  shown  in  in  table  X.  The 
budget  is  $74825.00  with  twenty  rive  stages  and  costs 
ranging  from  $10000.00  to  .50.  Due  to  the  length  of  the 
detailed  results  an  abbreviated  output  of  the  Marginal 
Analysis  solution  is  showL  in  table  XI.  The  return  is  2.853 
with  this  decision  set.  Note  in  the  last  line  of  the  table, 
the  split  in  the  bounds  indicates  the  heuristic  was  imple¬ 
mented  and  the  solution  may  not  be  optimal.  Examining  the 
abbreviated  output  of  the  Marginal  Analysis  program.  Table 
3.4  shows  that  the  budget  was  driven  negative  ten  iterations 
from  the  end.  Entering  the  same  parameters  and  budget  guan- 
tity  for  both  starting  conditions  of  the  Funnel  program 
produces  identical  results  shown  in  table  XII.  Comparing 
the  results  in  tables  XI  and  XII  it  may  be  seen  that  the 
return  and  decision  set  for  the  problem  are  the  same  for 
botn  programs.  Interestingly,  the  Marginal  Analysis  program 
was  able  to  correctly  select  ten  items  after  the  heuristic 
implementation  to  match  the  optimal  solution  produced  by  the 
Funnel  program.  The  heuristic  employed  by  Richards  and 
McMasters  appears  to  have  beer,  a  good  choice,  although  the 
solution  must  still  be  verified. 


TABLE  X 

Input  Parameters  for  a  Larger  Problem 


I 

DEMAND 

COST 

T 

ESSEN 

1 

2-000 

0.50 

1.90 

1.0 

2 

17.000 

3-00 

2.20 

3-0 

3 

0-250 

3.75 

2.50 

3.0 

4 

50-000 

5-00 

2.60 

1.0 

5 

1.000 

10.00 

1-60 

3.0 

6 

1-500 

12-00 

2.50 

1.0 

7 

18-000 

15.00 

1-50 

3-0 

8 

3-000 

20.00 

2-50 

1.0 

9 

0.250 

25.00 

1-50 

1.0 

10 

1.000 

25.00 

2.00 

1-0 

1 1 

2-000 

25.00 

1-50 

2.0 

12 

0-330 

50-00 

1.50 

3-0 

13 

0-100 

50.00 

2-00 

3.0 

14 

1-000 

61.00 

2-00 

2.0 

15 

7-000 

80.00 

2.30 

2.0 

16 

1-000 

91.00 

3-70 

3.0 

17 

2.500 

150.00 

1-75 

3-0 

18 

20.000 

210.00 

1-70 

3.0 

19 

3.000 

450.00 

2.70 

2-0 

20 

0.500 

500.00 

1.75 

1-0 

21 

1.  500 

710.00 

3.00 

2.0 

22 

6-000 

1000-00 

1-80 

1.0 

23 

2.000 

1500.00 

2.50 

3.0 

24 

12-000 

1500.00 

2.00 

2-0 

25 

0-500 

10000.00 

1.00 

1.0 

5  1 


Margina 1 

Analysis 

TABLE  XI 

Program  Results, 

- 

3udget=$74825. 00 

1  *  METHOD 

:  FIXED  EUDGET  = 

74825.00  * 

****************************** *********** ******  [ 

*********************************************,4  i  j 

STEP 

ITEM  SI 

RATIO 

BR  1  | 

1 

1 

1 

2-82236671 

74824.50  | 

2 

2 

1 

2. 14117622 

74821.50  '  1 

3 

2 

2 

2.08235264 

74818. 50  I 

4 

2 

3 

2.02352905 

74815.50  1 

5 

. 

. 

2 

4* 

• 

• 

4 

• 

4» 

1.9b470547 

• 

• 

74812.50 

440 

24 

21 

o!o004 187953 

23861 50 

441 

2 

55 

0.0004181664 

2383. 50 

442 

22 

9 

0.0003934940 

1383. 50 

443 

19 

1  1 

0.0003861645 

933.50 

444 

18 

40 

0.0003540027 

72J.50 

445 

4 

154 

0.0003520278 

7  18.50 

446 

24 

22 

0.0003425654 

-78  1.50 

446 

1 2 

3 

0.0003396191 

b68. 50  ■ 

447 

15 

23 

0.0003359483 

588. 50 

448 

7 

39 

0.0003278039 

573.50 

449 

8 

14 

0.0003014780 

553.50  • 

450 

20 

2 

G.  0002939 128 

53.50  J 

451 

4 

155 

0.0002806289 

48.50 

452 

2 

56 

0.0002605470 

45.50  i 

4  53 

10 

6 

0.0002369717 

20.50  ! 

454 

1 

12 

0.0002271836 

20.00 

455 

4 

156 

0.0002225503 

15.00  ! 

456 

7 

4G 

0.0002018004 

0.0  : 

i 

REALIZED  PERFORMANCE  =  2.8532 

i 

BCONDS 

< 

2.9796,  2.6762) 

i 

j 

i 

! 

Funnel  Program  Results,  Budget=$74825. 00 

THE  OBJECTIVE  IS  HSHT 
CONSTRAINT  ON  COST 

THE  PROBLEM  IS  TO  MINIMIZE  A  25  STAGE  PROCESS 
IN  IS  TC  BE  CHOOSEN  OPTIMALLY  BEISEEN  XN=7. 48250D+04 
AND  XN=7.48250D+04 

OPTIMAL  XN=7.43250D+04  OPTIMAL  BET0BN=2. 8 53 15D* 00 


N 

IN 

SN 

IN-  1 

25 

7. 4325D+04 

0.  0 

7- 4825D+04 

24 

7. 4825D+04 

2- 10D+01 

4.3325D+04 

23 

4. 3325D  *04 

6. 00D+C0 

3. 4325D+04 

22 

3.4325D+04 

9-OOD+OO 

2-5325D+04 

21 

2- 532SD+04 

6- OOD+QO 

2.  1065D  +  04 

20 

2. 1065D+04 

2. 00D+C0 

2.0065C+04 

19 

2-00650+04 

1. 10D+01 

1-51 15D  +  04 

18 

1. 51  15D  +  04 

4. 00D+01 

6-7  150D  +  03 

17 

6-7150D+03 

3. OOD+QO 

5.515GD+03 

16 

5- 51 50D+03 

8. 00D+00 

4.7870D+03 

15 

4-7370D+03 

2.30D+01 

2- 947CD+03 

14 

2- 947  0D  +  03 

5. OOD+CO 

2- 6420D  +  0  3 

13 

2.6420D+03 

2- OOD+OO 

2. 5420D+03 

12 

2. 5420D+03 

3- OOD+CO 

2-3920D+03 

1 1 

2.3920D+03 

7. OOD+OO 

2-2170D+03 

10 

2- 2170D+C3 

6. OOD+CO 

2.0670D+03 

9 

2- 0670D+03 

2. OOD+OO 

2-017CD+03 

8 

2.0170D+03 

1- 40D+01 

1-7370D+03 

7 

1-7370D+03 

4- OOD+Q1 

1.  1370C  +  03 

6 

1. 1370D+03 

9- OOD+OO 

1.0290C+03 

5 

1- 0290D+03 

6- OOD+OO 

9.690QD+02 

4 

9- 6900D+02 

1-56D  +  02 

1.8900D+02 

3 

1 - 8900D+02 

4- OOD+OO 

1. 7  400D  +  0  2 

2 

1.7400D+02 

5.60D+01 

6.0000D+00 

1 

6.0000D+00 

1-20D+01 

0.0 

-V  — J  ^  -V  — — T 


V.  CONCLUSIONS 

This  chapter  will  discuss  the  advantages  and  disadvan¬ 
tages  of  the  two  primary  solution  techniques  examined  in 
this  paper.  Various  starting  modes  of  the  Funnel  program 
will  be  looked  at  in  an  effort  to  develop  an  efficient 
approach  to  this  type  of  provisioning  problem. 

The  Marginal  analysis  program  provides  a  quick,  simple 
and  effective  approach  for  generating  a  solution  to  the 
previously  defined  provisioning  problem.  However,  the  solu¬ 
tion  is  not  guaranteed  to  b e  optimal  because  the  program  may 
implement  the  Richards  heuristic  to  consume  all  available 
resources.  The  Funnel  program  provides  the  optimal  solution 
to  the  same  problem  though  it  takes  more  computer  time  to 
operate.  The  DP  variant  has  the  capability  to  start  the 
problem  from  two  different  conditions,  1)  no  prior  knowledge 
and  2)  prior  knowledge  or  an  estimate  of  the  solution.  For 
Case  2,  utilizing  the  ENTER  function  of  the  Funnel  program 
allows  the  rapid  construction  of  the  input  data  file  used  by 
DP6  to  examine  the  Xi*  values  and  manipulate  the  window- 
boundaries  as  necessary.  If  the  entered  solution  set  is  not 
optimal,  the  Funnel  program  will  drive  the  decisions  to  the 
optimal  decision  set  in  subsequent  iterations.  Using  the 
timing  function  mentioned  in  Chap.  3  a  rough  approximation 
of  run  time  is  available  for  each  iteration  of  the  program. 
For  Case  1,  no  prior  knowledge,  the  total  estimated  time  for 
the  $74825.00  problem  is  nearly  3000  seconds.  For  Case  2, 
prior  knowledge,  a  near  optimal  starting  solution  for  the 
same  problem  is  driven  to  the  optimal  solution  in  approxi¬ 
mately  300  seconds.  This  901?  savings  on  computation  time 
makes  it  clear  that  an  approximation  method  as  an  input  to 
the  Funnel  program  would  be  a  better  approach  to  this  type 
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of  problem.  It.  should  be  noted  that  the  run  time  for  the 
case  2  start  is  related  to  the  accuracy  of  the  entered  deci¬ 
sion  set.  In  this  instance  the  entered  decision  set  was 
specifically  constructed  to  force  the  Funnel  program  to 
adjust  the  window  boundaries  in  all  stages  at  least  once.  A 
tetter  solution  than  the  one  entered  could  produce  a  much 
shorter  run  time. 

There  are  tradeoffs  in  the  Funnel  Program  for  computa¬ 
tional  efficiency.  Variations  of  the  multiplier  that 
creates  the  window  can  drive  the  individual  program  run  time 
directly.  A  window  that  is  too  large  wastes  computation 
time  on  resource  values  that  are  not  needed.  If  the  window 
width  is  too  small,  multiple  iterations  of  the  program  may 
be  required  to  compensate  for  the  increased  number  of 
boundary  violations  in  the  effort  to  reallocate  the  resource 
consumed  at  the  various  stages.  For  the  verification  of  a 
decision  set  provided  by  the  Marginal  Analysis  program  the 
multiplier  could  be  set  to  1.0  providing  a  very  narrow 
window  of  resource  to  be  examined.  Due  to  the  narrow  width 
of  the  window,  if  the  Xi*  value  changes  in  any  stage  the 
boundary  will  be  recorded  and  displayed  to  the  user  as 
violated.  This  notification  would  indicate  that  the  entered 
solution  is  not  optimal  for  the  problem. 

The  computer  time  expended  for  the  verification  alterna¬ 
tive  would  be  quite  low  even  for  large  budget  problems.  In 
the  25-stage  problem,  for  example,  approximately  one  second 
is  required  to  verify  the  correctness  of  the  optimal  solu¬ 
tion  if  it  is  entered  initially.  If  the  verification  alter¬ 
native  indicates  the  decision  is  not  optimal,  the  time 
required  to  drive  the  entered  solution  to  the  optimal  solu¬ 
tion  is  controlled  by  the  magnitude  of  the  changes  in  the 
decision  set  and  the  number  of  iterations  required  to 
re-allocate  the  resources  throughout  all  the  stages.  The 
user  must  make  a  value  judgement  on  the  importance  of  the 
optimal  solution  versus  the  additional  computer  costs. 
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In  conclusion,  the  emphasis  of  this  paper  is  to  create 
an  algorithm  to  determine  the  optimal  solution  in  a  quick 
and  easy  manner.  It  is  evident  that  a  combination  of  the 
Marginal  Analysis  and  Funnel  programs  is  an  efficient 
approach  for  solving  this  type  of  problem.  Using  the 
Marginal  Analysis  program  to  generate  a  best  estimate  of  the 
solution,  the  ENTER  f unct ion . witn  the  Funnel  program  can 
then  he  applied  to  confirm  the  optimal  solution  or  to  create 
an  input  data  file  that  DP6  may  use  to  find  the  optimal 
solution.  The  provisioning  problems  faced  by  NAVSUP,  as 
described  in  Chap  1,  may  be  solved  with  this  method  if  they 
can  be  formatted  in  the  recursive  DP  construct.  The  program 
DP6  in  particular  solves  only  the  minimization  of  MSRT  with 
a  constraint  on  cost.  Variations  on  this  program  coulu  be 
devised  to  solve  additional  objectives  like  supply  material 
availability (SMA ) ,  time  weighted  units  short(TKUS),  etc.,  by 
changing  the  computation  procedures  in  the  subroutines.  The 
constraint  however,  is  currently  restricted  to  the  linear 
form, 

22  Wi*Si<=E  i=1,...,N  (5.1) 

where  Wi  is  resources  required  per  unit.  Si  is  the  number  of 
units  of  item  i,  and  B  is  the  total  resource  available. 
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APPENDIX  A 
DP6  CODE 


This  appendix  contains  the  computer  listing  of  the  main 
program  and  subroutines  for  the  funneling  procedure.  Some 
of  the  subroutines  contain  code  for  objectives  and 
constraints  other  than  those  discussed  in  this  thesis. 
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